Dynamics of disordered heavy Fermion systems 
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Dynamics of the disordered heavy Fermion model of Dobrosavljevic et al. are calculated using an 
expression for the spectral function of the Anderson model which is consistent with quantum Monte 
Carlo results. We compute for three distributions of Kondo scales including the distribution 

of Bernal et al. for UCus-xPd^. The corresponding low temperature optical conductivity shows a 
low-frequency pseudogap, a negative optical mass enhancement, and a linear in frequency transport 
scattering rate, consistent with results in Yi-^U^Pda and UCus-^Pd^. 
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Introduction. Over the past four decades, the Fermi- 
Liquid paradigm has been the key to our understanding 
of metallic behavior. In Fermi-liquid theory we assume 
a 1:1 correspondence between the low- lying eigenstates 
of the interacting system to those of the noninteracting 
electron gas. This leads to a magnetic susceptibility that 
is weakly temperature dependent, a specific heat linear 
in T, and a low temperature resistivity that is quadratic. 

However, recently there has been a great deal of exper- 
imental interest in disordered heavy- Fermion compounds 
@ (e.g.Yi-zUzPds, UC^-^Pdz), where strong elec- 
tronic correlations preclude Fermi-liquid behavior. The 
non Fermi-liquid behavior in these compounds is charac- 
terized by a linear resistivity at low T, a logarithmic low 
temperature divergence of the susceptibility and the spe- 
cific heat coefficient, an optical conductivity with a low 
frequency pseudogap and a linear transport scattering 
rate at low frequencies |l] Several models have been 
proposed to explain these experimental results. Among 
them are theories based upon proximity to a zero tem- 
perature quantum critical point Q and those which ex- 
plain the impurity to lattice crossover effects in the multi- 
channel Kondo model || . 

A few years ago, Dobrosavljevic et al. investigated 
a system of dilute magnetic impurities in a disordered 
metal. Since disorder can give rise to a distribution of 
the local density of states of the conduction electrons, a 
distribution of Kondo scales P{Tk) is induced that could 
be singular enough to produce x ancl 7 that diverge as 
T — > 0, a strongly non-Fermi-liquid behavior. More re- 
cently, Bernal et al. fft have shown that this P{Tk) is 
not appropriate for the t/-doped heavy-Fermion systems. 
They propose an alternative spread of Kondo scales, and 
use it to calculate the thermodynamics of these alloys. 
Miranda, et al. || show that this distribution is consis- 
tent with a linear in T low temperature resistivity. 

Formalism. In this paper, we stay within this phe- 
nomenological framework and calculate dynamical quan- 
tities such as the optical conductivity and self-energy. 
We begin by concentrating on the Kondo regime of the 
Anderson impurity model, where d-electrons occupy the 



conduction band, and f-electrons provide the magnetic 
impurities for the spin-spin scattering processes. It is 
known that in the low temperature limit this model cor- 
responds precisely to the Fermi-liquid picture of Lan- 
dau, containing quasi-particle excitations from a ground 
state with relatively weak inter quasi-particle interactions 
(~ ksTjc). The electronic density of states has a reso- 
nance at the Fermi level, giving significant impurity con- 
tributions to the specific heat and magnetic susceptibil- 
ity. Following Doniach and Sunjic ||, Frota and Oliveira 
]n| argued that the Doniach-Sunjic form, modified to ac- 
count for the 7r/2 phase shift, should describe the shape of 
Kondo resonance. Their expression is in agreement with 
their results for the Kondo resonance obtained from nu- 
merical renormalization group calculations p"l| ], as well 
as low-temperature quantum Monte Carlo (QMC) re- 
sults analytically continued with the Maximum Entropy 
method (MEM) Jl|||. 

We can generalize their expression to finite tempera- 
tures and get 
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where A is the f-d hybridization energy, and Tk = 
(it/2) 2 Tk is the half-width of this resonance, Tk be- 
ing the Kondo temperature. We add a temperature 
dependent width 7 to the original expression of Frota 
and Oliveira p0[ , with 7 determined by fitting to QMC- 
MEM results pl2|. As shown in Fig.[l], this form contin- 
ues to exhibit remarkable agreement with the shape of 
the Kondo resonance (the low-frequency peak) obtained 
from QMC-MEM, even at finite temperatures. By com- 
parison with a wide range of Anderson impurity spectra, 
we were able to obtain the universal function j(T /Tk) 
for over three decades of T/T K - For T/T K < 0.3, we 
use ~f(T) — 4.52T, a value derived from Nozieres' Jji[ ] 
phenomenological Fermi-liquid description of the Kondo 
problem at low temperatures. The result for j(T ' /T K ) 
will be presented in a table of an upcoming publication. 

The Hilbert transform of Af((J, Tk) then gives the av- 
erage impurity t-matrix 
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t f (z) = JdT K P(T K ) J duo 



V 2 A f (uo,T K ) 
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where V is the f-d hybridization. Following Miranda, et 
al. H we use a dynamical mean-field approximation [ jl5| , 
which becomes exact in the limit of infinite dimensions, 
to calculate the lattice self energy from a concentration 
x of substitutional Kondo impurities 



E(w) 
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where £/(w) describes the average effective medium of the 
impurity. It is related to the average local host greens 
function G 



G(uo) 



through the relation 
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where JV(e) 



" e2 /** and we set t* = 10, 000K to 



establish a unit of energy and temperature. The solutions 
of Eqns. 1-5 then give the full lattice self energy. 

The knowledge of this self energy enables one to cal- 
culate physical quantities Jl6[ ] like transport coefficients 
and the optical conductivity. In this paper we concen- 
trate on the optical conductivity cr(w). It is measured 
in units of Co = e 2 ir/2fra, which with h/e 2 ~ 2.6 • 10 4 £1, 
varies between 10 _3 ...10~ 2 [(/zQcm) _1 ], depending on the 
lattice constant a. 

Results We choose three different distributions 
P(Tk), corresponding to strong ||, weak Q and a phe- 
nomenological disorder motivated by Miranda et al.'s ||] 
argument that the experimental P(Tk) should be rela- 
tively constant at low temperatures. Also consistent with 
Miranda we consider the distribution of Kondo scales as 
arising from a distribution of couplings between the con- 
duction and the f-electron spins P(J). For the strongly 
disordered sample, P(Tk) has the form M 



Psd(T K ) = (4tt) 



-1/2. 



T K ln(t*/T K ) 



exp{-0.251n 2 [0.217e _1 ln(t7T K )]} , (6) 

where we have used a bulk Kondo temperature T K = 
10 2 K (in most U-based HF systems T K = 10 2 K varies 
between 100-200 K). The weak disorder is characterized 
by a Gaussian distribution P(J) of width 2u — 0.01, 
where u is a disorder parameter j(| and a higher value of 
u corresponds to more disorder. This leads to 

P wd {T K ) 



v/OOTtt 0.217T K ln 2 (T K /t* 
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For the phenomenological spread of Kondo scales, we as- 
sume the form 



P P h(T K ) 



0.01 
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This distribution is not based on microscopies; it sim- 
ply satisfies the experimental criterion of constancy at 
low Tk and looks qualitatively similar to the P(Tk) for 
UCus-zPd,, g. Unlike the l/{T K \nT K ) divergence of 
the strongly disordered case, P p h{Ti<) has a finite num- 
ber of spins with Tk — 0. 

The nature of the disorder gives rise to different physics 
in each case, as is manifested in the functional form of 
ImE(cj) [Fig. 2] (these plots are for T = 0). For a weakly 
disordered system, ImSjui) has the form — c+uj 2 as uo — > 
where c is a constant |19|. This suggests a finite lifetime 
for the electrons at the Fermi energy at zero tempera- 
ture which makes it different from a normal (pure metal) 
Fermi liquid. But since there are no unquenched spins 
at T = 0, the system does form a local Fermi liquid, 
with a resistivity p(T) ~ p(0) — AT 2 . For the case of 
strong disorder, P S (i(Tk) is divergent at low Tk- Even 
though we are on the metallic side of the metal-insulator 
transition, a large number of spins with very low Kondo 
temperatures gives a non Fermi-liquid ground state. For 
this scenario, we find that IraS(o;) ~ — c + uu 1 / 4 at low u>. 

P p )i(Tk), the distribution of Kondo scales that is rel- 
evant to UCu5_ x Pda;, is intermediate between these two 
cases. It gives us an ImE(w) that behaves like — c + \u>\ 
as u> — > 0. We think that this behavior of ImS(w) 
can be understood through the following simple argu- 
ment. Oliveira's expression for the f-electron spectral 
function at T = can be expanded near uo — to 
give a form A/(a;,Tx) ~ 1 — a(u> /Tk) 2 , where a is a 
constant. Given that P p ii{Tk) ~ constant at low Tk, 
J Q Af(uj,Ti() P v }i{Tk) dTK gets its dominant contribu- 
tion from the region where Tk > uo. If we change the 
lower limit of the integral from to uo and make use of 
the fact that P p h{Tk) has a finite upper cutoff (~ 100K), 
ImE(cj) turns out to be — c + \uo\. In calculating the con- 
ductivity, the energy uo of the electron is averaged over 
a region of width kgT near the Fermi surface. This re- 
places |a>| by T, giving a resistivity linear in temperature 
pj which is observed experimentally in Yi-^U^Pda and 
UCu3.5Pdi.5 JlJj.This result is confirmed by direct calcu- 
lation of the resistivity (not shown). 

Fig. [3] shows the temperature dependence of the op- 
tical conductivity of the phenomenologically disordered 
system. Consistent with what is seen in Yi-^U^Pda and 
UCu3.5Pdi.5 0, the Drude-like peak is recovered when 
T > 100K, since spins with essentially all possible val- 
ues of Tk are participating in the dynamics at this tem- 
perature. The inset lists the zero temperature optical 
conductivities for the three distribution of Kondo scales. 
All are characterized by a vanishing Drude weight at low 
T, along with a finite frequency peak. As the disorder 
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is increased, this peak moves towards lower frequencies, 
concomitant with the decrease in the average value of Tk 
(which is different than T^). The Drude peak at T = is 
recovered from the weak disorder case if u — > in P(J), 
giving P wd {T K ) cx &{T K -T° K ) (Tg. is the bulk value), 
which takes us to the single Kondo scale physics. 

The optical conductivity of metals, even non Fermi 
liquid metals, is usually analyzed by rewriting it in a 
generalized Drude form || 



o-(w) 



1 



(9) 



47rr(u) + A(w)) ' 

where a(ui) = CTi(cj)+icr 2 (a') 0. We calculate the trans- 
port relaxation rate T(ui) and the optical mass enhance- 
ment (1 + A(cj)) for the three distributions of Kondo 
scales. However, as shown in Fig. 4(a), only the phe- 
nomenological distribution results in a linear in frequency 
zero temperature T(u), consistent with what is seen in 
Yi-zUzPds and UCu3.5Pdi.s |!J. In each case, we find 
that the low frequency optical mass enhancement 1 + A(0) 
is also negative (however, in the case of weak disorder a 
positive mass is recovered as u — > 0). 

r(aj) for the phenomenological distribution at several 
different temperatures is plotted in Fig. 4(b). Here, con- 
sistent with Degiorgi, et al. [0 we fit the relaxation rate to 
r(w) = r (1 - (T/T ) n - (w/wof) in the low frequency 
region. We find that T » 85 A' and lu w 0.09**, and 
are roughly constant in temperature. At high tempera- 
tures, r(w) and 1 + X(u>) (Fig. 4(c)) are weakly frequency 
dependent, 1 + > and n = 2.00, as expected for 
a Fermi liquid. Thus, as seen in Fig. 3, a Drude peak 
is recovered in <j\(w). As the temperature is lowered, 
r(w) and 1 + A(cl>) become strongly frequency dependent, 
1 + A(0) < 0, and n ~ 1, features which we believe should 
be viewed as characteristic of a non- Fermi liquid. A very 
similar sequence of features are seen in Yi-^U^Pds and 

UCu3. 5 Pdi. B 1- 

Within the formalism developed thus far, we can also 
calculate the magnetic susceptibility x(T) by disorder 
averaging Krishnamurthy's [fTl| universal susceptibility 
for a single Kondo scale. It is fairly straight-forward to 
compute the spin-relaxation rate in NMR (1/Ti) as well. 
These quantities will be addressed in an upcoming pub- 
lication. 

Conclusion. Within the Kondo disorder model , 
we calculate dynamics for some U-based heavy Fermion 
systems. We observe a linear resistivity at low T con- 
sistent with Ref. [8], the lack of a Drude peak and a 
low-frequency pseudogap in the real part of the optical 
conductivity, a negative low temperature optical mass, 
and a linear in frequency optical dynamical scattering 
rate. All these features are observed in Yi-^U^Pda 
and UCu3.5Pdi.5 Jl|,||. Thus, we conclude that the phe- 
nomenological distribution of Kondo scales model is suf- 
ficient to describe the dynamics of these disordered sys- 
tems. It is important to stress that Kondo disorder is not 



the sole possible explanation of non Fermi-liquid behav- 
ior in these systems [H||. In fact, it has recently been 
shown that the two-channel Kondo lattice model displays 
remarkably similar optical properties |l8[j . However, it 
remains to be seen if an appropriate two-channel Kondo 
model can accurately describe the transport and optical 
properties of dilute systems such as Yi_ x U x Pd3. 
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FIG. 2. Imaginary part of the conduction electron self en- 
ergy when T = and x — 0.2. The three different curves 
denote different distributions of disorder. The lowest one cor- 
responds to very weak disorder (P(Tk) — > as Tk —* 0) and 
at low uj has the form ImE(oj) cx — c + ui 2 , giving a local Fermi 
liquid with p(T) = p(0) — AT 2 . The highest one is for strong 
disorder (still metallic), with a u; 1 / 4 dependence as uj — > 0. 
The one in the middle corresponds to a phenomenological dis- 
tribution of Kondo scales suitable for the heavy-fermion sys- 
tems. It approaches the form ImS(w) cx — c + \u>\ as uj —* 0. 
This strongly hints towards a linear resistivity at low T in 
these compounds. 




FIG. 1. Spectral function for the f-electron in the Kondo 
limit of the Anderson model. The long-dashed line is the 
Doniach-Sunjic result jl^] (true for T = 0). The dotted lines 
represent the extension of the Doniach-Sunjic result to finite 
T which is used for calculating the dynamics in this paper; 
the solid lines represent QMC-MEM results. T/T K = 0.2, 
0.8, 3.2 and 12.8 for the curves from top to bottom. For 
T/T K < 0.3, Af(u,T h ) has a width 7 = 4.52T which comes 
from Fermi-liquid theory @. For higher T/T K , we adjust 7 
to fit the QMC-MEM data. 



-0.02 



-0.04 



-0.06 



-0.08 







\ 1 

\ \ 


/ / 
' / 
/ / 

/ / 

/ / strong 




/ . phenom 




/ / weak 



-0.2 



-0.1 




co/t 



0.1 




co/t 

FIG. 3. Optical conductivity for the phenomenological dis- 
tribution of Kondo scales at a few temperatures when x — 0.2. 
At very low temperatures a finite number of unquenched spins 
preempt the formation of a Fermi-liquid. The interesting fea- 
ture is the development of a Drude peak as we go from tem- 
peratures much below the bulk Kondo value (T K w lOOJf) to 
those much above it. The inset shows the conductivity for the 
three different P(Tk) when T = 0. The absence of a Drude 
peak is conspicuous in all three cases. 
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FIG. 4. (a) Frequency dependence of the scattering relax- 
ation rate T(ui) at T = for the three different P(Tk) when 
x = 0.2. Note that T{ui) is linear in uj, and therefore con- 
sistent with experiment ||, only for the phenomenological 
distribution of Kondo scales. T(uj) corresponding to weak 
disorder has an uj 2 behavior as uj — > 0, suggesting the for- 
mation of a local Fermi liquid, (b) T(uj) for P p h(Tk) at dif- 
ferent temperatures; in each case the solid line is a fit to 
the form F(uj) = To (1 - (T/T ) n - (uj/ujq)")- We see that 
up until 10K, we have a scattering rate that is roughly lin- 
early decreasing in uj and T. (c) Optical mass enhancement 
for P p ii{Tk) at different temperatures (symbols are the same 
as (b)). At low temperatures 1 + A(0) < indicative of a 
non- Fermi liquid. 



